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I. INTRODUCTION 

In the preceding publication [1] we started our analysis of the interactions between dark 
solitons of the parametrically driven nonlinear Schrodinger (NLS) equation: 

id T ^ + ld 2 x ^ + ^ -\^\ 2 ^ = h^* -i^. (1) 

This equation arises in a wide variety of physical contexts; see [l| for references. In Eq.flT]), 
h is the strength of the parametric driving and 7 is the damping coefficient. In the nondissi- 
pative limit, when 7 = 0, the equation has two coexisting stable soliton solutions, the Bloch 
and the Neel wall. In we considered forces existing between two Bloch walls and be- 
tween two Neel walls. The present work completes the analysis by classifying the Neel-Bloch 
interactions. The understanding of this nonsymmetric situation requires a mathematical for- 
malism different from the one used in jl|; this justifies the need for a separate treatment. 
Since the Bloch wall does not exist for 7 7^ 0, we consider here the nondissipative case only. 
This is another distinction from Ref. 

Our analysis of the interaction between Bloch and Neel walls will be based on the study of 
linearised perturbations about their (unstable) bound state. Mathematically, this analysis 
reduces to the construction of eigenf unctions of a Schrodinger-like operator, acting in the 
space of vector-functions, with the potential consisting of two well-separated nonidentical 
potential wells. When the two wells are infinitely far apart, there are four zero eigenvalues 
in its spectrum, with two associated eigenf unctions. As the wells (produced by the Bloch 
and Neel walls) are moved closer together, the degeneracy is partially lifted, with only 
two eigenvalues remaining at the origin. The question here is whether the two nonzero 
eigenvalues move onto the imaginary axis — which would be the case of stability, or onto 
the real axis — in which case the bound state is unstable. We use matched asymptotic 
expansions to show that the second is the case and construct eigenfunctions associated with 
the two real eigenvalues. 

The eigenfunction associated with the positive, unstable, eigenvalue contains the entire 
information on the character of interaction of two walls. We demonstrate that the outcome 
of the Neel-Bloch interaction depends on their left-right arrangement and the chirality of the 
Bloch wall. A Neel wall and a right-handed Bloch wall placed on its right will be moving in 
the same direction. If, however, we place a left-handed Bloch wall on the right of the Neel, 
the two walls will move in opposite directions — towards or away from each other, depending 



on the initial perturbation. The right-handed Bloch wall (on the left) and the Neel wall (on 
the right) will move in opposite directions, while a pair involving the left-handed Bloch on 
the left of the Neel, will move colinearly. 

After the eigenfunctions associated with two opposite real eigenvalues have been con- 
structed, the evolution of an arbitrary initial condition close to a pair of well separated 
Bloch and Neel walls is not difficult to predict. Treating this initial condition as a per- 
turbed Bloch-Neel bound state, its evolution will be determined by the projection of the 
perturbation on the bubble's unstable eigendirection. We illustrate this general approach 
by considering an example of initial condition in the form of a product of the Bloch and 
Neel wall. 

Usually one tries to understand the interaction of solitons as interaction of point-like 
particles; the particle description is physically appealing and mathematically lucid. We will 
show that a Bloch and a Neel wall can be treated as two classical particles. However the 
interaction between these two particles is anomalous in the sense that the Bloch wall being 
attracted to the Neel wall does not necessarily imply the reciprocal attraction of the Neel 
to the Bloch. This anomaly can be understood by considering the Neel wall as a particle 
with negative mass. The "wrong" mass sign arises very naturally if one recalls what the 
Neel wall really is: a localised depression, a patch of low density moving over a high-density 
background. The only reason why this property was not fully appreciated before is because 
earlier studies focussed on symmetric, Neel-Neel, interactions — which are, of course, non- 
anomalous. 

The outline of this paper is as follows. In section [III we introduce travelling Bloch and Neel 
walls, and describe the Bloch-Neel bound state. Section [TTT1 contains the main mathematical 
result of this paper, the asymptotic analysis of the splitting of the degenerate zero eigenvalue. 
In the next section (section HV|) we interpret the resulting eigenfunctions in terms of motions 
of the constituent walls. In section |V] we describe how the eigenfunctions can be used to 
classify the interaction of a pair of Bloch and Neel walls and apply this approach to a 
characteristic example. Finally, the main results are summarised in section IVII where we 
also interpret the interaction of the walls as interaction of opposite mass-sign particles. 
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II. MOVING BLOCH AND NEEL WALLS AND THE BLOCH-NEEL COMPLEX 

In this paper, we restrict ourselves to the undamped situation, 7 = 0. As in [1], we let 
V{X, T) = iAiP(x } t), x = AX, t = A 2 T } (2) 

with 



A = VTTh. (3) 

Equation (CD) becomes 

itj>t + ~ + fal> + j^V = 0. (4) 

This is the form of the parametrically driven NLS that will be used in this paper. The stable 
background solutions of Eq.(0| are V'flat = ±1- Without loss of generality we assume that 
h>0. 



The 



two topological solitons of (J4]) were introduced in Ref. [lj. One is the Neel wall 



iPn(x) = — tanh(x). (5) 

Note that we are introducing the Neel wall differently from [lj], with an extra negative sign 
in front of the tanh. This is done for later convenience. (In Ref.jl], we would refer to the 



solution (jSJ) as an antiwall.) 

The second topological soliton is the Bloch wall 



ip B (x) = tanh(Bx) ± iy/l - B 2 sech(Bx), (6) 

where 



A V 1 + h 

Equation ([6]) with the positive sign in front of the imaginary part describes the right-handed 
Bloch wall while in the case of the negative sign, the wall is said to be left-handed (see [q| 
for details). 

The Neel wall exists for all h > whereas the Bloch wall exists only for A 2 < |, i.e. 
for < h < |. Since we are interested in the Bloch- Neel interaction, the latter will be our 
region of consideration. Both the Bloch and Neel walls are stable in their entire regions of 
existence, < h < | and h > 0, respectively 9]. 



4 



The Bloch and Neel walls can be continued to nonzero velocity for all h values for which 



they exist 



9|. The moving walls of the form ip = ip(x — vt) are found as solutions of the 



ordinary differential equation 

- iwj> x + \p xx - HV> + ^ + JrV>* = 0. (7) 

The solution obtained by the continuation from the stationary Neel wall ip^{x), will be 
referred to as the "moving Neel wall" and denoted i/)n(v; x). In a similar way, the solution 
obtained by continuing ip B (x) in v, will be called the "moving Bloch wall"; to be denoted 
ip B (v\ x). The travelling walls are stable for all velocities v 

An important characteristic of solutions of Eq. (Tj0) is their field momentum: 

P = \ J (W x - W*)dx. (8) 

The momentum is conserved: dP/dt = 0. When we take the complex conjugate of ip(x,t), 
the associated momentum changes its sign: P[V>*] = —P[i(j]. Consequently, the momenta of 
two stationary Bloch walls with opposite chiralities are opposite. The right-handed Bloch 
wall has negative momentum, while the left-handed wall's momentum is positive. The 
momentum of the stationary Neel wall is, naturally, equal to zero. Figure [1] shows the 
momenta of stationary and travelling Bloch and Neel walls. These will be denoted Pb and 
P/v, respectively: 

P B = P B (v)=P[ip B (v;x)l (9a) 
Pn = Pn(v)=P[*I> n (v;x)]. (9b) 

According to Fig{TJ when the right-handed Bloch wall is continued into the region v > 0, it 
transforms into the moving Neel wall. In a similar way, when we path-follow the left-handed 
Bloch wall to large negative velocities, the corresponding branch turns into the branch of 
travelling Neel walls. Thus, the classification of moving solutions into Bloch and Neel walls 
is only sensible for sufficiently small v; for higher velocities, there is no qualitative difference 
between the two types of walls. 

One more observation with regard to Figfj] concerns particle properties of the two walls. 
While the stationary Bloch wall has positive mass, m B = (dP B /dv) v= o > 0, the mass of the 
Neel wall is negative: = (dP^ / 'dv) v= o < 0. This property will be crucial for the particle 
interpretation of the Neel-Bloch interaction (section IVl]) . 



p 


Bloch 
(left-handed) 


Neel 

^^"^ Bloch 
^ (right-handed) 


V 



FIG. 1: The bifurcation diagram for the stationary and travelling Bloch and Neel walls with h < | 
(adapted from Ref.[9]). The momentum ([8]) is used as the bifurcation measure. Solid curve indicates 
stable and dashed curve unstable solutions. Note that the mass of the Neel wall (defined as dP/dv) is 
opposite to the mass of the Bloch wall. In this plot, h = 0.05. 

In addition to the Bloch and Neel domain walls, Eq . (|1D possesses nontopological solitons. 
One such solution, arising for h = jz, is well known 
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12|: 



ib — 1 — -sech 2 (—} ± — tanh (—} sech ( 



(10) 



Recently, it has become clear that this solution is a member of a one-parameter family of 
solutions which exist for all h < | [13]. For each h, this family has the following analytical 
expression: 

4>b 



p + tag 



V 



where 



p — l _ e 4>i+W _ g2<fe+2/3 + e 4>l+2<f> 2 ^ 

q = 2(1 + B)e' h {l -e 01 ), 

and 

V — 1 + e^ 14 " 2 ' 3 + e 2</>2+2/3 + e 4 ' 1+24 ' 2 . 
In Eq.(TTTi). a is a sign factor: a = ±1; the exponents <fi 1 , 2 are 

(b 1 = 2(x + s), (b 2 = B(x-s), 

and (3 is defined by 

B = tanh (3. 
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(11a) 

(lib) 
(11c) 

(lid) 



The solution (fTTj) describes a bound state, or a complex, of a Bloch and a Neel wal 



with 



the parameter s characterising the distance between the centres of the two walls [13] . The 
choice s > corresponds to the Bloch wall on the right of the Neel wall; for s < 0, the Bloch 
wall is on the left. The solution ffTUj) corresponds to s = (and h = 1/15). The sign factor 
a determines the chirality of the Bloch wall bound in the complex: a = 1 implies that the 
Bloch wall is left-handed, while a = — 1 identifies a complex with a right-handed Bloch wall. 

Eq. ffTTj) can also be interpreted as a bubble of one phase (ip = —1) embedded in a 
background of a different phase (ip = 1); we will frequently be referring to this solution as 
a "bubble" Q. 

The analysis of the bound state (ITT|) will provide the understanding of the Bloch-Neel 
interaction. 

III. SPLITTING OF THE DEGENERATE ZERO MODES 
A. Stability problem for the complex 

To study the stability of the bubble (TTTj) . as well as to explore the phase space in the neigh- 
bourhood of this solution, we linearise Eq.fTjTJ) in the small perturbation 5ijj(x,t). Assuming 
the time dependence of the form Si[)(x,t) = [u(x) + iw(x)]e xt , results in the eigenvalue 
problem 

H(p = \J0, (12) 
where (p is a 2-vector made of the real and imaginary parts of the perturbation: 




J is a skew-symmetric matrix 



J= I ^ * h (13) 



and 7i is a hermitian operator: 




2KI 

U 2 + 3J 2 - ^# 



(14) 
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In the last equation, I stands for the 2x2 identity matrix, and 7Z and X are the real and 
imaginary parts of the solution ffTTT) : ipt, = TZ(x) + il(x). In what follows, we restrict our 
attention to the bound state with the Neel wall on the left of the Bloch wall [s > in Eq. fllip ]. 
The results for s < will be recovered by exploiting the symmetry ip(— x; — s) = ifj*(x; s) 
of the solution fill I) . Since the dynamics described by equation (J3J) are invariant under the 
reflection x — > —x, the evolution of a bubble with s = sq < will follow the same pattern 
as the evolution of the conjugate bubble with positive s = —s . 

For all values of s, the continuous spectrum of the operator 7i lies on the imaginary axis, 
with |ImA| > B, and does not give rise to instabilities. For s = oo, the separation of the 
Bloch and Neel wall in the bubble is infinite and so we have essentially two independent 
stationary (but potentially mobile) walls, each having two zero eigenvalues in its linearised 
spectrum. One of these stems from the translation invariance while the other one is associ- 
ated with velocity boosts of the corresponding wall. For finite s, only two of the four zero 
eigenvalues remain in the spectrum: one pertaining to the translation invariance of the com- 
plex as a whole and the other one associated with variations of the interwall separation. In 
this section we compute, perturbatively, the arising nonzero eigenvalues (and hence classify 
the stability of the bubble). We also construct the eigenfunctions associated with the real 
eigenvalues — these will provide insight into the evolution of the unstable bound state and 
nearby initial conditions. 

The lifting of the degeneracy of a repeated eigenvalue of the scalar Schrodinger operator 
with the potential comprising two identical potential wells with large separation, is discussed 
in the classical textbook 15]. This analysis is not helpful in our case, unfortunately, for 
three reasons: (i) our Ti operates on vector, not scalar, functions; (ii) the potential wells 
formed by the Bloch and Neel walls are not identical; (iii) the analysis in 15] postulates a 
particular form of the wavefunction on symmetry grounds, rather than deriving it within 
some perturbation formalism — as a result, the generalisation to the vector nonsymmetric 
case is not straightforward. 

Our treatment will be based on expanding the eigenfunction in the asymptotic series near 
the cores of the two walls and matching the resulting expansions in the overlap region x ~ 0. 
This approach builds on the asymptotic procedure used for the study of the stability of the 



travelling dark soliton 



16|- 



S 



B. Left expansion. 



We will do all our calculations for the case s > (i.e. for the Bloch wall on the right of 
the Neel wall). Furthermore, we will restrict ourselves to the right-handed Bloch wall only 
(a = —1). The other three possible combinations of a and Sign(s) will be commented upon 
at the end of section IIHI 

Let, first, x G (—00, Xq), where x > is a fixed value independent of s (so that x < s). 
In this region the real and imaginary parts of the bubble solution (fTTj) can be written as 



sinh[2(£-/?)] 
cosh 2 £ 



K = -tanh£ + \ 9 ^ c + 

J= 2(tanh£-5)e B(§ - 2s - /3) + (15) 

where 

£ = X + S + P, (16) 

and we have dropped terms of order e 3B ^~ 2s_ ^ and smaller. Equations ffl5l) can be seen 
as expansions of two functions of £, defined for —00 < £ < 00, in powers of e = e~ 2Bs . 
Accordingly, the operator 7i in Eq. (114ft expands as 

H = H N + eH i i ) +e 2 H% ) + ... (17) 

Here, 7i N is the unperturbed linearised operator of the Neel wall centred at £ = 0, i.e. 
Eq. (jl4p with 1Z = — tanh£ and 1 = 0. Guided by the results of (our own) numerical 
analysis, we assume that A is of order e: 

A = \ e. 

(As we will see, this asumption leads to a self-consistent perturbation scheme.) This implies 
that the eigenfunction ip in the eigenvalue problem (fl2|) can also be expanded in powers of 
e: 

0(0 = aiJ' N (0 + e0i(O + e 2 M0 + (18) 

Here ip' N = d^ipN, and ipN = (Hn,Zn) = (— tanh£, 0). The coefficient a is arbitrary at this 
stage [a = 0(1)]. 

Substituting (jTTT)- (TL8T) into (fl2l) and equating coefficients of like powers of e, yields, at 
order e 1 : 

H N {p\ + aH { ^if' N = a\ Jif' N . (19) 



To solve (119]) . we note that ip N is a generalised eigenvector associated with the zero eigen- 
value: 

n N if N = -Jtf' N . (20) 

Here ipx is considered as a solution of the form ip = ip(x — vt) of equation (JTj), while the 
overdot indicates differentiation with respect to velocity v (and not time): ip N = ( d v ipN ) 
(This will be our convention until the end of section ITO1 ) Also, we know that ip' b = d x ipb 
and d s ^b are zero modes of the perturbed operator (Tbfl) where ipb(x) = (7^>2T) is the bubble 
solution (ITTj) . and — oo < x < oo. Using (JT5J) , we write 



where 



= $v + £ (xx - ft) + (x 2 - m) + (21) 



X 1 (0 = 2e^ ) sech 2 e( ° ], (22) 
y x (0 = 2Be B ^\t^i - B) ( ° ] , (23) 

^.^--E^a/i,, (24) 



cosh f I o 



(25) 



Substituting expansions (TTTT) and (l2"Tj) into 7Y^£ = H S V>& = 0, the order gives a useful 
identity: 

H N x 1 + H ( i ) ft N = ] (26) 

we also note that 

H N y! = 0. (27) 
Using (1201 and (T2"6"j) . we can solve Eq. tlTD"]) in the class of functions bounded as |£| — >■ oo: 



<Pi(0 = -a^N + a>X 



i- 



The emerging vector-function = aip' N +e0i + ... decays, exponentially, both as £ — > — oo 
and £ — » +oo and hence it is intuitively clear that it cannot describe the behaviour of the 
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eigenfunction of the bound state in the region under consideration (x < 0). Indeed, the 
eigenfunction should include a term growing as £ — > +00 and representing the tail of the 
Bloch wall situated in the region x > 0. Therefore we need to add to (p a solution of the 
equation TCnV — which decays as £ — > —00 but grows as £ — ► 00. This solution is already 
available [see (1271) ]; it is given just by Eq. (l2"3"|) . Note that the other two linearly independent 
solutions of equation TCnV — (other than ip' N and yi), are growing as £ — > —00; these are 
clearly not acceptable for our purposes. 

Finally, the full order-e 1 perturbation is —a\ip N + aexi + Ci^i, and the eigenfunction (p 
becomes, in the region x < Xq ( "left region" ) : 

<f(C) = a$ N ~ a\i/J N + aexi + Cm + 0(e 2 ), (28) 

with 2/1 as in fl23|) and the constant C\ to be found from the matching condition at a later 
stage. (Here we are implicitly assuming that this constant will be of order e. If, however, it 
is of order e 2 or higher, the term Cy\ will only appear at higher orders of the expansion.) 



C. Right expansion. 

Let, now x G (— Xq,oo). (We remind that x > is a fixed value independent of s; 
x < s.) Here, the real and imaginary parts of the bubble are given by 



U = tanh (Brj - 0) - e~ 2 ^ + 

cosh 2 (Br) - (5) 



x = sech fi 2 coshjBr]) c _ 2(r?+2s) , . 

cosh(Br7 - 0) cosh 2 [Br] - (3) '"' 1 J 



where 



T] = x — s (30) 

and we have dropped terms of order e~^ v+2s \ Eqs. (1291) can be seen as expansions of the 
functions TZ(rj) and I(f]) (with — 00 < rj < 00) in powers of \x = e~ 4s . Accordingly, the 
operator TC expands as 

n = n B + + ^nf + .... (31) 

Here, Hb is the unperturbed linearised operator of the Bloch wall centred at r] — (3/B, 
Eq. (fl4j) with TZb — tanh (Br/ — (3) and X B = sech (3 sech (Brj — (3). Note that the expansion 
parameter /i is, in general, incommensurate with e and hence it is not a priori obvious what 
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the expansion of the eigenfunction <p{rj) will be. Letting <p{rj) = ip' B + <p, where (p has some 
order of smallness, and substituting into (fT2l) . we get 



H B 0=\ o ejft B , (32) 

where we have dropped terms of order /i because /i is smaller than e (and even e 2 ). From 

fl32|) it is clear that the leading-order correction to the translation mode of the free-standing 

Bloch wall is of order e (and not \i or /i 1 / 2 as one might have been tempted to think). 

Recalling that ^>b(v) = \9 v iPb) is a generalised eigenvector associated with the zero 

V / v=0 

eigenvalue of the free-standing Bloch wall: 

= -J4>'b, 

the localised solution of Eq. (l32l is given by <p = —\ip Bl and so (the localised part of) the 
eigenfunction (p can be written as 

0{r ] )=H' B -\bif B + O{e\ii), (33) 

where b is an arbitrary constant of order 1. 

The eigenfunction (|33l) decays, exponentially, as rj — ► — oo and therefore, cannot represent 
the (x > 0)-behaviour of the eigenfunction of the bound state. (For there is no connection 
of such an eigenfunction to the region x < where the Neel wall is located.) In order to 
obtain the correct behaviour, we need to add a solution of equation Tisz = which decays 
as r] — » +oo but grows as r\ — > — oo. The equation has only two solutions that decay as 
rj — > oo; one is the translation mode 

fa = CJZ' B ,1' B ) = (b sech 2 (B V -0),-B sechp ~ | ^ , (34) 

and the other one can be found by expanding the exact null eigenfunctions of the bubble, 
ip' b = d x ipb and d s ijjb, in powers of \x (in the region x > —xo). We have 

d s j b = -fa + fi (z x + ^i) + 0(/i 2 ), (35) 



where 



m - T s smHWv) I (36) 

cosh 2 (^-/5) \ 2cosh(^) ! 
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and 

SM- 1 ?! <«**)■ w 

Substituting (l3Tj) and (1351) into 7^0j, = 'Hd s i ) b = 0, we obtain T^b^l = which means that 
Eq. (l36l) gives exactly the solution we are looking for. 

Finally, the correct behaviour of the eigenfunction <p in the region x > — x (the "right" 
region) is given by 

0(r l )=b^ B ~\b^ B + Dz 1 + O(e 2 ), (38) 
where z\ is given by Eq.( l36i) and the constant D is to be fixed later. 



D. Asymptotic matching. 



Equations (12811 and ( 1381) give the leading terms in the asymptotic expansion of the eigen- 
function in the regions x < Xq and x > —Xq, respectively. Before proceeding to the next 
order of the expansion (which will give us the eigenvalue A), we need to make sure that the 
two expansions match in the overlap region — xq < x < Xq. 

In what follows, we will need the asymptotic behaviours of the generalised eigenvectors 
ip N and ij) B . The solution ip N of the equation (1201) can be constructed in quadratures. The 
|£| — > oo asymptotic behaviour is straightforward from the quadrature: 




Mt) - ~P( B + !) e IC ' I I . ( 39 ) 
where the factor 

1 POO 1 JD Icy 

» = m=m L^ {B ' tanhfleB ^ = i-b>m*b/2) > a (40) 

The rj — > ±oo asymptotics of the solution ip B of equation Hb^b = ~J4 ) 'b are gi ven by 

u ( u ± e-\ B ^\ \ 

Ml) - 2|B al . 41 

I W±e -2|B„-/?| J 

where the constants w± and w± are easily found by direct substitution: 

45(1 + B)e~P 8B 
U± = T 4-B* ' W± = 4^- 
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Using (1391) . the £ — > +oo-asymptotic behaviour of the expansion (1281) is given by 

_ 4ae -2(x+ S +/3) \ / 

y y p(5 + i) e -B(*+»+/J) 

+ 8a £ e^) ( ° ] +d ( ° | . (12) 

where we have substituted x + s + (3 for each £. The rj — > — oo- asymptotic behaviour of the 
expansion (13"Hj) is 



(B-2)(a:-s)-2/3 



4j g e 2B(x- S )-2/3 \ / u _ e B(x-s)-f3 \ / _4 e -2(*-»)-2/9 

- Afc +^ I - (43) 



5e 



Here we have replaced each r] with x — s. 

Equating coefficients of the exponential e Bx in the bottom rows of fj42|) and fH3l) deter- 
mines the constant G\\ 

d = be. (44) 

Equating coefficients of the exponential e~ 2x in the top rows of (142]) and (j43l) and the 
exponential e^ -2 ^ in the bottom rows, fixes the constant D: 

D = apt. (45) 

(Both top and bottom rows lead to equivalent equations, so only one parameter is fixed.) 
Since D turns out to be smaller than A and even A 2 , we can drop the term Dz\ from (p\ and 

fa- 
There is a term in ( 1421) which does not have a matching partner in (j43l) . and the other way 
round, there are terms in (1431) which do not have counterparts in (1421) . Consider, first, the 
term proportional to e~ Bx in Eq.( T42l) . This exponential does not have a partner in Eq.(f4"3l); 
however we can add a matching term —a\tp N (£) with £ = t] + 2s + (3 to the expansion (l38l) . 
This term will be of order e~ 4Ss for 77 ~ and hence will arise only at the next, e 2 — , order of 
the expansion. In a similar way, the exponentials e Bx in the top row and e 2Bx in the bottom 
row of (1431) do not have counterparts in the expansion ( 1421) . This can be taken care of by 
adding the term —b\ip B {rj) with 77 = £ — 2s — (3 to Eq.( f28l) . The top and bottom rows of 
this term will appear only at the order e 2 and e 3 , respectively. Finally, the last unmatched 
exponential e 2Bx in the top row of (1431) will acquire a matching partner if we add e 2 y2 to 
Eq.fTSHl). where j7 2 (0 is given by Eq. (r25|) . 
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E. Secular equation. 



To identify the constants a, b and the eigenvalue A we proceed to the next order of the 
perturbation expansion. As before, we treat the regions x < and x > separately. For 
— oo < x < x , the order e 2 gives 

H N (p2 + Hx<pi + aHni)' N = \qJ<Pi, (46) 
where <pi is as in Eq. (1281) . The solvability condition of Eq. fHBI) is 

(j' N , wgVi) + a (jJ' N , H^n) = A (A, M) > (47) 
where (, ) stands for the L 2 -scalar product of two vector-functions of £: 

(fj)= r (/v + /V)de. 

J —oo 

To obtain (j47|) . we have made use of the identity (iI>' n ,Hn<P2) — {.02i'Hn'>P'n)i an d the fact 
that H-niP'n = 0- However, since ^ should include terms which have exponential growth 
as £ — > +oo, the validity of this identity (and hence of the equation (I4T1) ) may be under 
suspicion. To reassure that the solvability condition is indeed correct, we note that the 
second-order perturbation should have the form 

eV2(0 = ~b\J B + e 2 y 2 + C 2 y x + e 2 2l (48) 

where C 2 is a constant of order e 2 , and the function (p 2 is bounded as |£| — > oo. Note that 
the growing terms in the right-hand side of fT4Hl) (the first, second and third terms) grow no 
faster than e 2B ^ whereas ip' N decays as e -2 ^ when £ — > +oo. Therefore, when doing each of 
the integrals (-0^, Hn^Pb)^ C0jy ^NVi) an d C0jv> 'HnVi) by parts, the boundary terms vanish. 
Consequently, we have (^' n ,Hn02) = ( ( P2,'Hn'^n)'i Q-E.D. 

Eq. (|47p can be simplified considerably if we use the identities arising at the order e 2 of 
equations Hip' = and Tlips = with Ti. expanded as in ( flTl) and t/>', as in ( J2TT) : 

^%+^!; ) xi+^ivx2=o, 

W$j7l + = o. 

Taking the scalar product with ip' N , these become 

fa,H$tf N ) + (f N ,H$&) = 0, 

(A,<Vi)=0. (49) 
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Using Eqs.(j26D, (HI) and (ggjl. Eq.gTD simplifies to 

aA 2 J^) - teA = 0. (50) 



Noting that ( tp' N , Jip N j = (1/2)P N and evaluating the integral yp' N , Jyij , Eq. floTJ]) becomes, 
finally, 

yJVa + 2e\pB 2 (l - B 2 ) e~ Bp b = 0. (51) 



Here Pn is the momentum of the Neel wall, Eq. (19bl) . 

Turning to the region — xq < x < oo, the eigenfunction has the expansion 

0(r ] )=b^ B -Xb^ B + e 2 2 + O(X 3 ^), (52) 

where (p 2 consists of a bounded part (p 2 an d a part that grows as r\ — > — oo: 

eV 2 fa) =A 2 £ ? -Aa^ v (£). (53) 

The order A 2 of the expansion of H</? = A </</? gives 

(aV 2 - aA^) = -X 2 bJ$ B . (54) 

The solvability condition is 

- aA = -A 2 6 J^ B ) , (55) 

where, this time, 

At first glance, the scalar product on the left-hand side of (1551) is zero as TCb^Pb = 
0. However, ip N grows as r\ — > — oo: -^tv ~ e _B,? , while does not decay fast enough 
- in fact, it decays at exactly the same rate: $ B ~ e Br) . Therefore, {^p' B , r H,B'^N S j 7^ 
ip N , TCb^'b ) ■ (Note the difference from our analysis of the neighbourhood of the Neel wall 



where ip' N decayed faster than the exponentially growing terms in <p 2 would grow, and hence 
(^P' n ,Hn02J = (^^jT^n^'n) was indeed true.) Making use of the asymptotic expression 
fl39|) and taking care of the boundary term, we get, instead: 

($ B , H B j N ) = 2epB 2 {\ - B 2 )e- Bl3 . (56) 
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Finally, Eq. (l55|) becomes 



2e\pB 2 {l - B 2 )e~ BI3 a + —P B b = 



A 2 



(57) 



A 



0. 



where Pb is the momentum of the Bloch wall, Eq. (l9a|) . 

Eqs. (l5Tl) and ( 1571) constitute a system of two linear equations for two unknowns, a and 
b. The associated characteristic equation is fourth-order in A: 

XP N AepB 2 (l- B 2 )e~ Bf3 

AspB 2 {\- B 2 )e- Bf3 -XP B 

The equation has two zero roots corresponding to the translation symmetry of the bound 
state as a whole and variations of the interwall separation; there is also a pair of real roots 
of opposite sign: 

B 3 e -Bf} 



A = ±2tt — c e - 2g *. (58) 

(Note that P N and Pb are opposite in sign — see Fig{TJ hence we have a positive quantity 
under the square root). This is the first result of our analysis. The second result is the 
relation between a and b, the coefficients of the leading terms in the expansion of (p: 



a 
b 




(59) 



Here the plus corresponds to the positive eigenvalue and minus to the negative one. (That 
is, the "unstable" eigenfunction has a and b of the same sign, both positive or negative, 
whereas the "stable" eigenfunction has the coefficients of the opposite sign.) In what follows 
it will be convenient to choose the normalisation such that a = \[p~b > for either sign of 
A, while b = \/—Pn > for A > and b = —^/—Pn < for A < 0. Thus, the normalised 
eigenfunctions are: 



, A , \JP B (ip' N - \X\^ N + exi) + y-PNsyi + 0{e 2 ), 
0s = \[Pb fe + |A|Av + exi) - V-iWi + 0{e 2 ] 



(60a) 



for x < 0, and 



<Pu^\-Pn[^b-Wb)+0(e 2 



P N ($ B + \\\j B )+0(e 



(60b) 
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for x > 0. 

We conclude our calculation by commenting on the case of the complex involving a left- 
handed Bloch wall (cr = +1) and on the situation where the Bloch wall (right- or left-handed) 
is on the left of the Neel wall (i.e. s < 0). These cases do not require a separate treatment. 
Indeed, let (u(x),w(x)) be the eigenvector associated with the eigenvalue A (positive or 
negative) of the bubble (TTTT) with s > and o = — 1. Changing o — > —ex changes the sign 
of the off-diagonal entries of the matrix (fT4l) ; hence the vector (u(x), —w(x)) will give the 
eigenfunction associated with the eigenvalue —A for the bubble with a = 1. Therefore, the 
"unstable" eigenvector for the bubble with s > and a = +1, obtains from the "stable" 
eigenvector of the bubble with s > 0, a = —1, just by changing the sign of its bottom 
component. To find the eigenfunctions for the Bloch wall on the left of the Neel, we note a 
symmetry of the bubble solution ffTTj) : ipb{x, —s) = ip%(—x, s). According to this symmetry, 
changing s — > — s in Eq.f llll) is equivalent to replacing x — > —x and taking complex conjugate. 
Consequently, the vector (u(— x), — w(— x)) will serve as the eigenfunction associated with 
the eigenvalue —A for the bubble with a = — 1 and s < 0. Finally, (u(— x), w(— x)) and A 
are the eigenvector and the associated eigenvalue pertaining to the situation with a = +1, 
s < 0. 



F. Numerical verification 

The above perturbation results were verified numerically. (We used Fourier expansions 
of (p over 120 positive and 120 negative harmonics, on the interval [—40, 40].) In agreement 
with our expectations, for s = oo we have found four zero eigenvalues. Two of these remain 
at zero as s is decreased from infinity. The other two zero eigenvalues move to the real 
axis; one becomes positive and the other negative (FigfSJ). Their magnitudes are equal. 

, with H 



131 . The 



This is a general property of linearizations of Hamiltonian systems of the form ( TT2 
hermitian and J as in f|T3l) ; it follows from the conservation of symplectic areas 
eigenvalues return to zero as s — > 0. These are the only eigenvalues with nonzero real part 
in the spectrum of 7i; hence for any finite, nonzero s, there is exactly one unstable mode. In 
addition, there are two pure imaginary eigenvalues detaching from the continuous spectrum 
as s decreases from infinity. These do not reach the origin; they remain on the imaginary 
axis and hence do not cause instability. 
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FIG. 2: Two discrete real eigenvalues of the eigenvalue problem (TL2!) . The solid lines depict eigenvalues 
found numerically while the dashed line gives the perturbation approximation (l58l) . This plot is for the 
case where the Bloch wall bound in the complex is right-handed; changing the chirality of the Bloch 
wall switches the two branches, but leaves the overall shape the same. This figure pertains to h = 0.2; 
for other h, the functions A(s) look similar. 

In addition to the numerical solutions of the eigenvalue problem ( ITSl) . we performed 
numerical simulations of the full time-dependent nonlinear Schrodinger equation (j3]). Simu- 
lations were carried out using a split-step pseudospectral method under periodic boundary 
conditions (see e.g. jisl)- Typically, we used an interval (— L/2, L/2) with L = 60, although 
for small h and large s (i.e. when the walls bound in the bubble are far apart and decay 
slowly to the background), L = 120 was necessary. The timestep was set at At = 2.5 x 10~ 3 . 
The code is stable for At < 7r _1 (L/A^) 2 , so for L = 120 we were able to use N = 1024 
modes, while for L = 60, we were limited to N = 512. The spatial resolution Ax = 2nL/N 
was the same in both cases. Results of numerical simulations are presented in the next two 
sections, concurrently with predictions of the asymptotic analysis. 

IV. INTERPRETATION OF THE EIGENFUNCTION 

The eigenfunction associated with the eigenvalue A > admits interpretation in terms of 
the motions of the constituent walls. Let, for defmiteness, s > and a = — 1, and note that 
the dissociating bound state can be described by 

f(x, t) = ip N (v N {t); x + s + 5s N {t) - x {0) ) + ... (61a) 
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FIG. 3: The eigenfunctions of the bubble involving the right-handed [a = — 1, (a,b)] and left-handed 
[a = +1, (c,d)] Bloch wall. The eigenfunction in (a,c) is associated with the positive and that in (b,d) 
with the negative eigenvalue. In this figure, s = 5 and h = 0.2. The real part is shown by the solid 
line, while the imaginary part is dashed. 



in the region x < 0, and 

ip(x, t) = ip B (v B (t); x — s — 8s B (t) - x (0) ) + ... 



(61b) 



for x > 0. Here ip^iy^x) and ip B {v]x) are stationary solutions of Eq.([7]) with small v 
obtained by continuation from 4>n(x) and ^b(x), respectively, and includes the part 
of the perturbation which cannot be reduced to the translation and velocity boost of the 
corresponding wall. Note that since we have not specified this part, the parameters vn, 
vb, 8sn, and 5sb are not defined uniquely in (|6T|) ; to fix these parameters we need to 
restrict the "..."-part in some way. To do this, we note that Eq. fl6TT) can be represented as 
i^bix) + Sijj(x,t), where ipb{x) = is the stationary bubble ffTTl) with some s and x^°\ 

and a small perturbation 5i()(x,t) can be written in the form 



5ip(x, t) = v N (t)ip N + 5s N (t)ip' N + z7(x, t), 
5ip(x,t) = v B {t)ipB ~ Ss B (t)ip' B + u(x,t) 
in x < and x > 0, respectively. Here 

ip N = d v i/; N (v;x + s - x (0) ) ; 

v=0 

^ N = dJ N (0;x + s-x (0) ); 
i) B = d v ip B (v;x - s - x {0) ) ; 

v=0 

ip' B = d x ip B (0; x-s- s (0) ). 



(62a) 
(62b) 



20 



(We remind that the overdot stands for the partial derivative w.r.t. v, not t, here.) For the 
given s and x^°\ we can fix vn, vb, 5sn, and Ssb by requiring that the "remainder" u(x,t) 
be J-orthogonal to the subspace spanned by the velocity boosts and translations: 

^ rO 

(u, Jip N ) dx — (z7, Ji/)' N ) dx = 0, (63a) 

-oo J — oo 

/■OO . />oo 

/ (u, J$ B ) dx= (u, Jf B ) dx = 0. (63b) 
Jo Jo 

Here (, ) denotes the scalar product in the two-dimensional Euclidean space: (/, g) = f\gi + 
]i92- From (1631) we have 

vn = — — 

(i/; N , Ji)' N ) 

and similar expressions for vb, 8sn, and 5s B . One more consequence of the constraints (|63|) 
is that v(x, t) is linearly independent from ip N and ip' N in the region x < and from ip B 
and ip' B in the region x > 0. (Indeed, if we assumed z7 = Ci^yv + C 2"*/V an d substituted 
this in H63aj) . then using the fact that (i[) N , Ji[) N ) = \Pn ^ 0, we would immediately get 
ci = c 2 = 0.) 

The perturbation ( 1621) can be expanded over solutions of the equation ( [121) : 

(ty = M u e xt u + .M s e~ At & 

+ M^' b + J\f d s if b + Q(+) e ^°V (+) + Q ( " ) e-^°V ( ~ ) 

/oo /*co 
Q( 1 )(A;)e^^ 1) cifc+ / Q {2) (k)e iuJt ^ ] dk. (64) 
-oo ./ — oo 

In Eq.( l64l) . (p u and (p s are the eigenfunctions associated with discrete real eigenvalues A 
and —A, respectively, where we set A > 0; ip' b and d s ipb are the two zero modes while 
and cj)^ are eigenfunctions associated with pure imaginary discrete eigenvalues ±icu . 
Finally, 0^ and are solutions of the continuous spectrum, = iujJ(f)^' 2 \ with 

uj 2 = \{k 2 + A)(k 2 + B 2 ), and cu(k) > for k > and w(Jfe) < for < 0. 

As t grows, the expansion (1641) tends to A4 u e xt (p u , where (p u {x) is given by ( |60|) . This 
should be identified with the large-t behaviour of (1621) . Using linear independence of vectors 
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i/j' n , ip N , V in x < and vectors ip' B , ip B , V in x > 0, we get 

v N (t) -> -XM u \fKe x \ 



v B (t) -> -\M U ^/ -P N e Xt ; (65) 
5sjv(t) -> M u \[p~ B e x \ 



5s B {t) ^ -M U \J -P N e xt (66) 

as t — > oo. Equations ( l65i) and ( l66i) are consistent in the sense that the velocity of the 
wall determined from the deformation of its shape coincides with the velocity defined by the 
position of its centre: 

Note that the derivative Pat is greater, in absolute value, than P B (see Fig{T]). Conse- 
quently, Eq. fl65l) implies that the Bloch wall arising from the dissociation of the bubble, 
always moves faster than the emerging Neel wall: \v B \ > \vn\- The same Eq. (1651) implies 
that the corresponding velocities and accelerations of the Bloch and Neel walls will be colin- 
ear. Consequently, the walls emerging from the decay of the bubble with s > and a = — 1 
will be moving in the same direction. (The direction will of course be determined by the 
initial perturbation.) This conclusion is in agreement with direct numerical simulations of 
Eq.®. (See Fig. 11(a)). 

It follows then from equation (1651) that if M. u > 0, the two walls will be moving to the 
left, with the Bloch wall catching up with the Neel wall. This was indeed seen in simulations 
(Fig. Ufa)). (In our numerical simulations, the choice Ai u > was accidental; the initial 
perturbation of the bubble was entirely due to the discretisation errors and hence beyond 
our control.) On the contrary, if Ai u < 0, the walls will be moving to the right, with the 
Neel wall lagging behind. 

Changing o — > —o swaps around the top components of the eigenfunction (p u associated 
with A > and eigenfunction (p s pertaining to the negative eigenvalue —A. Therefore, the 
"unstable" eigenfunction of the bubble with a = +1 involves coefficients a and b of the 
opposite signs. As a result, the fragments of its decay — the Neel and left-handed Bloch 
wall — will be moving in opposite directions. The simulations confirm this (FigJH(c)). 

As we have already mentioned, the choice of the coefficient A4 U was beyond our control 
in the simulations. For s = 5 and a — +1, our discretisation induced M. u > 0. However, for 
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FIG. 4: The development of the instability of the bubble (fTTi) . Shown are the curves of constant real 
part. The wide and narrow trails pertain to the Bloch and Neel walls, respectively. Panels (a) and (b) 
correspond to bound states involving right-handed Bloch walls (a = — 1), with s = 5 and s = —5, 
respectively. The evolution of the bubbles with the left-handed Bloch walls (a = +1) is depicted in 
panels (c) and (d), also with s = 5 and s = —5, respectively. No perturbations were added to the 
stationary initial condition (fTTj) "by hand"; the initial (uncontrollable) disturbance was merely due to 
the discretisation of the equation. The collision of two walls in (a) produces a breather. This breather 
moves so quickly that the contour plot seems to indicate the presence of two separate breathers whereas 
there is in fact only one. This effect is due to a sparse time-sampling of the profile which we had to 
resort to for better visualisation. The image in (a) has been generated by saving a profile once only 
every 200 time units. In these plots, h = 0.1. 

a slightly different value of s, viz. s = 4.9, the "unstable" eigenfunction (p u was seen to be 
excited with the coefficient M. u < (Figj5]). 

The effect of the change s — >• — s amounts to changing o — > —a and x — > —x; therefore, 
the dissociating bound state of the right-handed Bloch wall on the left of the Neel wall 
will produce walls moving in opposite directions. [Simulation shown in FigJU^b).] Finally, a 
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FIG. 5: Evolution of the bound state involving a left-handed Bloch wall (cr = +1), with s = 4.9. This 
initial condition is almost identical to the initial condition that led to Figfjjc); however, the direction 
of motion of the walls is now exactly opposite. The reason, of course, is that M. u was positive in 
the perturbation of Fig|4jc), whereas here, the "unstable" eigenfunction (p u is excited with a negative 
coefficient. As in Fig|4] in this plot h = 0.1. 

left-handed Bloch wall (a = +1) and a Neel wall placed on its right (i.e. s < 0) will move 
in the same direction. If (as it happened in our simulations) M. u is > 0, the Bloch will be 
leaving the Neel wall behind [Figj4](d)]. 

In the case of converging walls, the result of the collision in all cases is the formation 
of a spatially localised, temporally oscillating object (which we refer to as a breather), 
propagating over the constant background ip = 1. An asymptotic expression for the breather 



was derived in the previous submission [l]; here, we simply note that it persisted indefinitely 
after it was formed in our simulations. 



V. THE BLOCH-NEEL INTERACTION 



A. The method. 



Our knowledge of the phase space in the vicinity of the unstable bound state of the 
Bloch and Neel walls can be used to make conclusions on the evolution of particular initial 
configurations of the two walls. Any initial condition comprising well separated Bloch and 
Neel walls can be represented as 

ip(x, 0) = ijjb(x) + 5ijj(x, 0), 
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where if)},(x) is the bubble (TTTT) with some suitably chosen separation s, centred at some point 
a;*- -*, and 5ip(x,0) is a small perturbation. This initial perturbation can be expanded as in 
(I64p where we just need to set t = 0. In order to find the coefficient Ai u with which the 
unstable eigenvector (p u enters the perturbation, we note that the J-product of the vector 
(p s with any solution of Eq. (fT2l) except (p u is zero: 

(<p s , J(p s ) = (<p s , Jip' b ) = (<p s , Jd s ip b ) 
= (ft, J$&) = (ft, J$>) = (ft, J$f) = 0. 

[This is a simple consequence of the hermiticity of the operator 7i in Eq. fll2l) : if TiC,! = ^i^Ci 
and HC2 = A2JC2, then (A x + A 2 )(Ci,^C2) = 0. Hence (Ci, ^2) = unless Ai = -A 2 .] 
Therefore, 

M u = (67) 

((p s , Jifu) 

The sign of M. u determines the direction of colinear motion of the two walls. In fact, 
the denominator in (|67|) is positive and the sign of Ai u is determined just by the sign of 
(0s, Jdijj). Indeed, using the representation ( 160al) for (p s (x) } <p u (x) and invoking the identity 
( 1BTI]) . one can check that the integral over the negative semiaxis is zero to the leading order 
in e: 



(<p a , J0 u )dx = 0(e 2 ). 

-00 

We remind that (, ) is the scalar product in the two-dimensional Euclidean space: (f,g) = 
fi9i + /2<?2 whereas (, ) stands for the L 2 -scalar product: yf(x), g(x) J = J^ o {f,g}dx. On 
the other hand, using the expansion ( I60bl) for < x < 00, we obtain 



(tp 8 , J0 U ) dx = -XP N P B , 
which is positive for A > (see FigHj). Thus, (<p a , J0 U ) = -\P N P B + 0{e 2 ) > 0. Q.E.D. 



B. Example. 

As a characteristic example, we consider the initial condition of the form 

ip(x, 0) = -ip N (x + Xi)ip B (x - xi) (68) 

with some X\ > 0. In applications, this product function is often used as an approximation 
for two well-separated dark solitons, in the same way as two weakly overlapping bright 
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solitons are usually approximated by their sum. (See e.g. our previous submission [l| where 
this type of ansatz was employed for the variational analysis of the Neel-Neel and Bloch- 
Bloch interactions.) Assuming that the Bloch wall is right-handed and that x\ is large, we 
have 

Re^(x, 0) = -tanh£ + 2e 2B ^~ 2xi) tanh£ + 

ImV>(x,0) = 2sech ( 5e' B(5 - 2:Cl) tanh£ + ... (69a) 

in the region x < 0, and 

Re^(x, 0) = tanh(5r/ - /3) - 2 e -^+^-P/ B ) tanh(^ - /?) + 
kni/)(x,0) = sech/Jsech^ - 0) - 2 sech/3 e -2(i+2x 1 -^/B) sech ( B7] - p) + ... (69b) 

in the region x > 0. Here we have defined £ and 77 according to 

£ = x + xi, 77 = x-£i + ^. (70) 

£5 

Note that the leading, O(e ), terms in (1691) coincide, formally, with the leading terms in the 
asymptotic expansions of the bubble, Eqs. (jT5l) and (!29|) . However the definitions of £ and 
77 in flTDl will not, in general, be consistent with the definitions of £ and 77 for the bubble, 
i.e. Eqs.f ll6p and (1301) . To achieve the consistency, we introduce an additional translation 
parameter a;( ) in the definitions ([TBI) and O30p : 

£ = x + s + /3 } rj = x - x (®)- s . (71) 

Thus we allow for translations of the bubble to ensure the coincidence of the leading terms in 
(I6"9"|) and (fT5|) . (T2"9~j) . Equating (170]) to (I7T|) we obtain parameters of such a reference bubble: 

•--IH). ---IH)- 

Subtracting the bubble solution with parameters (I72p from the initial condition (1681) . we 
obtain the following expressions for Re Sip = Keip(x,0)—TZ(x) and Im^ = Im^(x, 0)— X{x): 

e 2(B-l)« 

Re5V(^) = £V 2B/3 sinh(2/3) 5— + 

cosh £ 

Im # (x) = 2eBe~ Bp — + . . . (73a) 

cosh 4 
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in the region x < 0, and 

2(B-1) V 

Re5ip(x) = fie- 213 sinh(2/3) = + 

V ; P V cosh 2 (Br) - (3) 

_ 2 cosh(Bri) - e" 2/3 sech/? cosh(Br] - (3) 

Imdw (x) = 2ue ' ^- h ... 

^ y 1 P cosh 2 (5r/ - /3) 

(73b) 

for x > 0. Here e = e~ 2Bs and /i = e _4s ; 7?. and X are the real and imaginary parts of the 
solution (fTTT) : ipb = TZ + il. 

In order to evaluate the integral (<p s , JSijj) in Eq. (l6T|) . we note that the positive semiaxis 
of x gives a contribution of the order e -( 2 + B ) s , This is exponentially smaller than the con- 
tribution of the negative semiaxis (which is of the order e~ 2Bs ) and hence can be neglected. 
Using (I60al) and (I73aj) . we get then 

(<p s , J5$) = 2e^J-P N Be- Bp ! e^-^sech 3 ^ + ... 

J — oo 

which is positive. Hence M. u > 0, and, according to Eqs. fl65|) - fl66|) . both walls will be moving 
to the left. 

The direct numerical simulations of Eq.(jl]) verify these conclusions. Fig J6] presents simula- 
tions of the initial condition with x\ = ±6 and h = 0.1, for both chiralities of the Bloch 
wall. (These results are representative of all h provided is sufficiently large). When the 
Bloch wall is right-handed and X\ is positive, i.e. when the Bloch wall is on the right of the 
Neel, both walls move to the left (Figj6](a)) — precisely as our analysis predicted. Using 
symmetries of the initial condition fl68|) and the evolution equation (j3J), one can readily check 
that the situations shown in Figjf3(b-d) are in agreement with the asymptotic analysis as 
well. 

VI. CONCLUDING REMARKS 
A. Energy considerations 

Equation (jlj) is conservative, with the energy integral given by 

e = \ J {foM 2 + H 4 - JrM 2 W + m 2 ] + l} dx. (74) 
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FIG. 6: Typical results of the evolution of an initial condition of the form (|68l) . Panels (a) and (b) 
correspond to the right- handed Bloch wall, to the right {x\ > 0) and the left {x\ < 0) of the Neel 
wall, respectively. Panels (c) and (d) correspond to the left-handed Bloch wall, to the right and the 
left of the Neel wall, respectively. The Bloch wall moves towards the Neel wall in all cases. The Neel 
wall moves towards or away from the Bloch wall depending on whether it is on the left or right of the 
Bloch wall, and on the Bloch wall's chirality. After the collisions in (b) and (c), a fast-moving breather 
is formed; however, for visual clarity, this solution is not shown. Here, h = 0.1. 

The bubbles, being time-independent solutions of Eq.(|4]), must render the functional (!74|) 
stationary: 5E = 0. Since they depend on s smoothly, it follows that dE/ds = 0: the 
energy of the bound state is independent of the distance between its constituents. This 
seems to suggest that the binding energy of the Bloch and Neel walls is zero, implying the 
non-interaction of the walls. 

However, our analysis has revealed that the walls do, in fact, interact. We have shown 
that all bound states with the parameter s ^ are exponentially unstable against the decay 
into constituent walls. Exponential growth of the intersoliton separation (see Eq. (1661) ) is a 
clear manifestation of a nonvanishing interaction between the walls — for noninteracting 
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walls, the separation would only grow linearly. As for the zero binding energy, it can be 
reconciled with the nonvanishing interaction by noticing that the bound Neel wall acquires a 
small static imaginary part (see Eq. (jl50 ). (The bound Bloch wall also acquires a stationary 
perturbation, but it is exponentially smaller than the perturbation of the Neel wall.) To 
leading order, the force between two walls is made up of the force between their real parts 
and the force between their imaginary parts while the force between the real part of one wall 
and imaginary part of the other one, is of the second order of smallness [3]. The imaginary 
part of any nontrivial stationary solution decays as 0(e~ Bx ), while the real part decays as 
0(e~ 2Bx ) or 0(e~ 2x ). Consequently, the force caused by the overlap of the imaginary parts 
of the soliton tails is stronger than the force caused by overlapping real parts. Thus although 
the amplitude of the imaginary excitation of the Neel wall is exponentially small for large 
s, the interaction of this excitation with (the imaginary part of) the Bloch wall is enough 
to balance the force between the Bloch wall and the "naked" Neel wall. If we decrease s, 
and the two walls in the complex are pulled closer, the increase of (the absolute value of) 
the binding energy of their real parts is offset by storing more energy of the opposite sign 
in the imaginary excitations; the total energy remains invariant. A similar mechanism was 
described by Ostrovskaya et al in the context of the dark-bright solitons of the undriven 
vector nonlinear Schrodinger equation [20] ; there, the imaginary excitations were referred to 
as "solitonic gluons" . 



B. Walls as particles 

Next, we need to explain, qualitatively, the anomalous behaviour apparent in Figs|4](a,d), 
and [6] (a,d), where the velocities and accelerations of the two interacting walls are seen to 
have the same (rather than the opposite) directions. The key observation here is that 
the Bloch and Neel walls, considered as point-like particles, have masses of opposite signs 
(see Fig{T]). One can easily conceive a simple model system of positive- and negative- mass 
particles, which exhibits the observed phenomenology. Let X\ and x 2 be the coordinates of 
the two particles, with z = X\ — x 2 , and let mi > and m 2 < be their masses. Consider 
the Hamiltonian 

*=2^'-« 2+ 2^-^ < 75 > 
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where pi,p2 are the momenta of the particles, and the functions fi = fi(z) and fi = /^(z) 
decay, exponentially, as \z\ — > oo. The equations of motion are 

oxi mi m 2 

dH ■ (7^ 

P2 = -77— = -Pi, (76b) 
ox 2 

£i = 77— = — (pi - fi), (76c) 
dpi mi 

an i 

X2 = 77- = — (P2-/ 2 ), (76d) 

where the overdot indicates differentiation with respect to time it) (and not v as earlier in 
the text). The prime stands for the derivative w.r.t. z. 

For any z = zq, the equations ( 176|) have a fixed point, describing an unstable bound state 
of two particles: 

Pi = fi(z ), P2 = f2(z ). (77) 

This fixed point is an analogue of the bound state of two walls, Eq.flTTJ). Linearising Eqs. flTB"]) 
about the equilibrium point (I77|) . and letting <5p 12 (t), Sxi^it) oc e xt , we get a pair of nonzero 
real eigenvalues 

A (±) =±^=[/{(.o) + /^(.o)]. (78) 

(There is also a pair of zero eigenvalues resulting from the overall translations xi^ — > 
%i,2 — and the freedom in choosing zq.) Working out the associated eigenvectors and 
substituting them into the linearised equations fl76cl) - (176dl) . we find 

mi5ii = ± — Sz, (79a) 

vM 

K (±) 

\m2\Sx2 = Sz, (79b) 
yml 

where are coefficients dependent on mi j2 and /{ 2 ( 2; o)- 

If f[(z ) + fzizo) > 0, the unstable eigenvalue is A^, and we keep the top sign in fl79al) . 
This is the case of anomalous behaviour: the velocities Sx\ and 8x2 as well as accelerations 
Sxi and S'x2 have the same sign and so the particles move in the same direction, like the 
walls in FigsJU (a,d), and [6] (a,d). If f[(zo) + /^(^o) < 0, the positive eigenvalue is A^ - -' 
and the "unstable" eigenvector is given by the bottom equation in (179aj) . In this case the 
two particles have opposite velocities and opposite accelerations; this corresponds to the 
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"normal" interaction of the Bloch and Neel wall seen in the panels (b,c) of Figs JH and [6j (It 
is fitting to note here that the anomalous behaviour does not mean that the two particles 
violate the third Newton's law. The force exerted by the particle 2 on particle 1, —dH/dx\, 
is exactly opposite to the back reaction force, —dH/dx2-) 

C. Conclusions 

Finally, we summarise the results of this investigation. 

The essence of our approach is to consider the interacting Bloch and Neel walls as a 
perturbation of their unstable stationary complex. The interaction between the two walls is 
characterised by the eigenfunction associated with the positive eigenvalue in the spectrum 
of the linearised operator (evaluated at this stationary solution). 

Using matched asymptotic expansions, we have evaluated the real eigenvalues for the 
Bloch-Neel complex and constructed the associated eigenf unctions. The structure of the 
"unstable" eigenfunction for the complex consisting of a right-handed Bloch wall on the 
right, and a Neel wall on the left, indicates that the walls emerging from the decay of this 
"bubble", will be moving colinearly, i.e. in the same direction. This rule determines the 
evolution of a general configuration of a right-handed Bloch wall on the right and a Neel wall 
on the left — as long as the walls are sufficiently far away from each other. Using symmetry 
properties of the Bloch-Neel complex, we can also predict the type of motion (colinear 
or antilinear) for other chiralities and mutual arrangements of the walls. The asymptotic 
analysis is in agreement with direct numerical simulations of the interacting Bloch and Neel 
walls. 

Although the unstable eigenfunction determines the type of motion of two interacting 
walls, knowing the structure of this eigenfunction is insufficient to know which direction 
this motion will take. (Depending on whether the unstable eigenfunction is excited with a 
positive or negative coefficient, the colinearly moving walls may travel to the left or to the 
right. Similarly, in those cases where the unstable eigenfunction sets the opposite direction 
of motion for the walls, they may travel either towards or away from each other.) The actual 
direction of the colinear or antilinear motion depends on the particular initial condition and 
can be determined by the projection of the corresponding perturbation of the stationary 
complex on its unstable eigendirection. We have evaluated this projection for an initial 
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condition in the form of a product of the Bloch- and Neel-wall solutions, and verified the 
conclusions of the asymptotic analysis numerically. 

Finally, we have interpreted the anomalous interaction of the Bloch and Neel walls as a 
dynamics of two interacting particles, one with positive and the other with negative mass. 
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